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Eigenvalue correlations in QCD with a chemical potential 
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We discuss a new Random Matrix Model for QCD with a chemical potential that is based on the symmetries 
of the Dirac operator and can be solved exactly for all eigenvalue correlations for any number of flavors. In the 
microscopic limit of small energy levels the results should be an accurate description of QCD. This new model 
can also be scaled so that all physical observables remain at their /i = values until a first order chiral restoration 
transition is reached. This gives a more realistic model for the QCD phase diagram than previous RMM. We also 
mention how the model might aid in determining the phase diagram of QCD from future numerical simulations. 



It is well known that at zero chemical potential 
the low energy modes of QCD can be described 
by a chiral effective Lagrangian 

C=^-Tv[d^Ud^} -|Tr[A/([/ + [/t)] . (1) 

In a finite volume and at low energy when 
1/TOtt 3> L (also known as the e-regime) the ki- 
netic term can be ignored and the theory reduces 
to simply an integration over the zero momentum 
part of the SU(iV/) matrix U that parameterizes 
the Goldstonc manifold. The results then do not 
depend on the dynamical details of QCD but rely 
solely on the symmetries. In this limit the chiral 
effective theory can also be expressed as a Ran- 
dom Matrix Theory (RMT) with the appropriate 
chiral symmetry |T|. The partition function is 
given in terms of an integral over all (N + v) x N 
complex matrices A 
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with a Dirac matrix 
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The choice of a complex matrix corresponds to 
a SU(N C > 3) gauge field with v the topological 
charge. In the microscopic limit of N going to 
infinity with Nm held fixed the chiral RMT is 
equivalent to the static part of the chiral effec- 
tive Lagrangian. The RMT form of the effective 
theory, however, is better suited for studying the 
eigenvalues of the Dirac operator. 



Here we want to examine the extension of the 
chiral RMT to include a baryon chemical poten- 
tial. In this case a Hermitian term /170 is added 
to the anti-Hermitian p. The eigenvalues of the 
Dirac operator become spread out in the complex 
plane and this gives rise to a sign problem that 
hinders numerical investigations. 

One model considered previously by Stephanov 
|5] used the Dirac matrix 



V = 





L4t + fj,l 



iA + fj,I 




(4) 



where the chemical potential 11 is multiplied by an 
identity matrix. This has the correct symmetries 
but the fixed form of the chemical potential term 
makes it very difficult to calculate the eigenvalues 
analytically. Also it has been shown that this 
model has an unphysical \i dependence when it is 
studied at larger li 

A new model for QCD with a chemical poten- 
tial was introduced recently [|]. This uses a dif- 
ferent form for the Dirac matrix 



V = 







iA + fiB 




(5) 



where the chemical potential part is now also 
modeled with a Gaussian random matrix B. This 
model still has the correct symmetries and can 
be solved for all eigenvalue correlations analyti- 
cally. Additionally we can remove the unphysical 
\x dependence at zero temperature giving a more 
physical model for chiral symmetry breaking with 
a chemical potential. 
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Figure 1. Example of a direct measurement of 
a first order chiral phase transition. To deter- 
mine [i c accurately one needs several data points 
around p c . 




Figure 2. By extrapolating the effective parame- 
ters a and p to zero one could determine fi c at a 
first order phase transition without direct simu- 
lations at larger chemical potentials. 



By choosing an appropriate parameterization 
for the matrices A and B and then integrating 
out all angular variables the new model can be 
written in terms of the complex eigenvalues z k of 
© as g] 

N f N N f 

Z oc J] m v j / |A(z 2 )| 2 11 w(z k )dz k ]J(z 2 + m?) 
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with the Vandermonde determinant 



A(z 2 ) = Jft, 



(7) 



i<3 



and a weight function containing a modified 
Bessel function 
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This representation makes it easy to solve for all 
eigenvalue correlations for any number of flavors 
using the orthogonal polynomial method. The 
results are given in 0] which we won't repro- 
duce here. In the microscopic limit the eigen- 
value correlations from the RMT are expected to 
agree with QCD due to universality. Recent lat- 
tice studies support universality in the eigenvalue 
density at small chemical potentials 



We now examine the behavior of this model at 
larger chemical potentials. If we look, for exam- 
ple, at the result for the partition function for one 
flavor given by a Laguerre polynomial 
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we see that the fi dependence comes in only 
through the combination cr/(l — /i 2 ). AtT = 
we expect no fi dependence of thermodynamic 
quantities until a first order phase transition is 
reached. We can remove the fj, dependence com- 
pletely by setting a(fi) — a (1 — p 2 ). We then 
find that all thermodynamic observables such as 
the chiral condensate in the chiral limit 
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(10) 



now do not depend on the chemical potential until 
a phase transition is reached. This therefore pro- 
vides a reasonable effective model for QCD with 
a chemical potential at zero temperature. 

We can also use this as an effective model at 
nonzero temperature. To do this we first include 
an explicit scale for /i in the Gaussian weight for 
the matrix B so that the matrix weights are 

exp(- cr iVTr [A^A] - P NTt [B f B] ) . (11) 

This gives the same results as replacing p with 
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P\J a j p in previous expressions. The chiral con- 
densate in the chiral limit is then 
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From this we see that there are two possible sce- 
narios for the chiral condensate to vanish as p is 
increased. One is that a goes to zero while p re- 
mains nonzero. This will be first order only if a 
has a discontinuous jump to zero, otherwise it is 
second order. Another possibility is that both a 
and p reach zero at the same point. In this case 
the transition can be first order even if a and p 
are continuous. 

The results of the effective model may be very 
useful to future numerical studies. For example 
in order to locate the possible first order chiral 
phase transition at large p with a direct simu- 
lation one would need several measurements of 
the condensate near the transition as illustrated 
in figure 2] This is after performing the neces- 
sary thermodynamic extrapolation to find a sharp 
transition. Due to the severe sign problem, per- 
forming the necessary simulations at such a large 
p is likely to remain prohibitive for a long time. 
An alternative approach is to measure the effec- 
tive parameters a and p at smaller p and then 
attempt to extrapolate them to zero. The critical 
point is given by when a reaches zero. If p also 
reaches zero at the same spot then it could be a 
first order transition. This allows one to find the 
critical point without going to very large p. Of 
course current simulations still have difficulty at 
the moderate p needed to extrapolate, so it will 
likely still be a while before the numerical simu- 
lations are feasible. 

The most general form of effective theory for 
chiral symmetry breaking with a baryon chemical 
potential can be obtained by letting <r(T, p) and 
p(T, p) be functions of T and p. The parameters 
a and p can be determined from numerical simu- 
lations by fitting to a quantity such as the eigen- 
value density or the valence quark mass depen- 
dence of the chiral condensate. The latter is likely 
to be a better observable since it is smoother and 
should be easier to obtain with higher statistics. 
Then by mapping out the parameters in the T-p 
plane one could end up with a phase diagram like 





2nd order 


o->0 






\ 1st order 




a jumps to \ 




or \ 




a --> \ 




p --> 1 



Figure 3. Possible phase diagram. 

the example in figure |31 

We now have a new Random Matrix Model 
for QCD with a chemical potential that can be 
solved for all eigenvalue correlations for any num- 
ber of flavors. In the microscopic limit this pro- 
vides exact expressions that should agree with fi- 
nite density QCD due to expected universality. 
This model can also be used as an effective theory 
for chiral symmetry breaking in the presence of a 
chemical potential. By considering the effective 
parameters in the model to be general functions 
of T and p this may provide a useful way to help 
map out QCD phase diagram. 
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